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j^ , Abstract 

We study the short-time dynamics of a mean-field model with non-conserved order parameter 
r^ • (Curie- Weiss with Glauber dynamics) by solving the associated Fokker-Planck equation. We obtain 

o: 

CJ . closed-form expressions for the first moments of the order parameter, near to both the critical and 

spinodal points, starting from different initial conditions. This allows us to confirm the validity of 
^sg ■ the short-time dynamical scaling hypothesis in both cases. Although the procedure is illustrated 

l/-^ . for a particular mean-field model, our results can be straightforwardly extended to generic models 

o, 

I I I with a single order parameter. 

o 
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I. INTRODUCTION 

Universal scaling behavior appears to be an ubiquitous property of critical dynamic sys- 
tems. While initially believed to hold only in the long time limit, it was realized during 
the last decade that the dynamical scaling hypothesis can be extended to the short-time 
limit [l|. This is accomplished by assuming that, close to the critical point, the nth moment 
of the order parameter obeys the homogeneity relation 

m(") (t, r, L, mo) = 6" "'^/"m^") {b-'t, h^'^r, L/b, 6^mo) , (1) 

where t is time, r is the reduced temperature r = {Tc — T)/Tc, L is the linear system size, 
mo is the initial value of the order parameter and 6 is a spatial rescaling parameter, /x is a 
universal exponent that describes the short-time behavior, while /3, u, and z are the usual 
critical exponents. When m-o << 1 we recover the usual dynamic scaling relation, from 
which a power law relaxation at the critical point (for instance, in the magnetization n = 1) 
m{t) ~ t^P/^^ when L ^ 1 and t ^ 1 follows. This is the critical slowing down. 

On the other hand, the short-time dynamics (STD) scaling properties of the system 
depend on the initial preparation, i.e., on the scaling field rriQ. Setting h = t^/^, from Eq. ([1]) 
one obtains for small (but non-null) values of t^^^rriQ in the large L limit 

m(t, r, mo) ~ mo t'F{e'''T), 9 = l^—lll, (2) 

z 

Hence, at the critical point r = an initial increase of the magnetization m ~ mot^ is 
observed. For the second moment n = 2 the dependency on mo can be neglected when 
mo ^ 1. Since m^'^^ ~ A^^^ = L^*^ in the large L limit (d is the spatial dimension), one 
obtains at the critical point 

(t) ~ A^-4'^/^-2^/"^ (3) 



m 



(2) 



The short-time universal scaling behavior has been verified in a large variety of critical 



systems both by renormalization group (RG) calculations 







and Monte Carlo (MC) 



numerical simulations J3|-l5|. The hypothesis also applies when the system starts in the 
completely ordered state, i.e., mo = 1. In this case it is assumed that the homogeneity 
relation 

m(") (t, r, L) = b-^'^/'^m^'^^ {b-H, b^'^r, L/b) (4) 



holds even for short (macroscopic) time scales. Hence, in the large L limit we have 

m{t) = t-^/'^"G(ti/^V), (5) 



and taking the derivative of logm, 



d\ogm(t, t) 



~ t^/"'. (6) 

r=0 



dr 

While the scaling hypothesis starting from the disordered state is supported both by 
numerical simulations and RG, its validity for an initial ordered state relies up to now only 
on numerical simulations. 

Recently, numerical simulations have shown that the short-time scaling hypothesis ([1]) 
holds not only close to a critical point, but also close to spinodal points in systems exhibiting 
a first-order phase transition, both for mean-field and short-range interactions models |6|. 
This is particularly interesting, because it suggests the existence of some kind of diverging 
correlation length associated to a spinodal point. Since the proper concept of spinodal 
in short-range interactions systems is still a matter of debate (see Ref. [61] and references 
therein), a deeper understanding of the microscopic mechanisms behind the observed short 
time scaling could shed some light on this problem. One way of achieving this goal is to 
look for exact solutions of particular models. A first step in that direction is to analyze 
mean- field (i.e., infinite-range interactions) models, for which the concept of spinodal is well 
defined 6|. That is the objective of the present work: we analyze the exact STD behavior 
of far from equilibrium mean- field systems with non-conserved order parameter. 

Non-equilibrium phenomena in physics and other fields are commonly studied through 
Fokker-Planck equations (FPEs). In particular, non-equilibrium dynamical aspects of phase 
transitions can be analyzed by means of the FPE associated to the master equation de- 
scribing the microscopic dynamics |7H9|. In fact, this tool was proved to be useful in the 
description of the relaxation of metastable states [7;], finite-size effects |10| or the impact 



of fluctuations in control parameters 



ll| . and have been considered for mean- field spin 



models [T], |l2| and coupled oscillators [105, amongst many others. 

As soon as the degrees of freedom of the system can be reduced to a few relevant ones, a 
low- dimensional FPE can be found. Although this description is suitable for properties that 
do not depend on the details of the dynamics, or for mean-field kinetics, many conclusions 
are expected to hold in more general instances. 



For a single order parameter m, the FPE for its probability P = P{m^ t\mQ, 0) is 

dtP = l-dmDiim) + dmmD2{m)] P = LFp{m)P , (7) 

where the drift and diffusion coefficients are determined by the Hamiltonian and the partic- 
ular dynamics (e.g., Glauber or Metropolis). 

Following this stochastic approach, here we study the scaling of the short-time relax- 
ational dynamics in the vicinity of critical and spinodal points. In first approximation, 
the drift Di{m) (= —dV/dm) is generically linear in the vicinity of a critical point and 
quadratic in the spinodal, following the quadratic and cubic behavior of the drift potential 
V, respectively. Meanwhile, typically in various models, the noise intensity D2{m) scales as 
e ~ 1/A^ 3, |lO|. Therefore, although we will present the STD for a particular spin model, 
our results can be straightforwardly extended to more general mean-field ones. 

II. FORMAL FPE SOLUTION AND MOMENT EXPANSIONS 

The formal solution of the FPE ([71), for the initial condition P{m, 0|mo, 0) = 5(m — mo), 

is y 

P{m, t\mo, 0) = e*^^^("')5(m - mo) . 

The average of an arbitrary quantity Q{m) can be derived directly from the FPE, by 
means of the adjoint Fokker-Plank operator Lpp{'m) = Didm + D2dmmi as follows 



{Q){mo,t) = I Q{m)P{m,t\mo,0)dm= Q{m)e^^^''^"'^6{m - mo)dm 
S{m - mo)e*^^^("')g(m)t/m = e*^^^("«)g(mo) = 
= ^[4p(mo)]'=Q(mo)tV^!- (8) 

fc>0 

Therefore, the first two moments of the order parameter are 

(m) = mo + D,t + ^[DiD[ + D2D';]t^ + . . . , 
(m^) = (m)2 + 2D2t + [2D2D[ + D^D'^ + D2D!^]f + ... , (9) 

where Di, D2 and their derivatives are evaluated in mo- Notice that if Di and D2 are not 
state-dependent, the expansion up to first order is exact. 



Alternatively, evolution equations for moments can be obtained by integration of Eq. ([7]), 
after multiplying each member of the equation by the quantity to be averaged, that is 

^^^ = n {m''-^Di{m)) + n{n - 1) (m^-^/^alm)) . (10) 

For n = 1 we have 

d(m) , , X, , . 

-^ = {D,{m)). (11) 

Eqs. (fTO|) lead in general to a hierarchy of coupled equations for the moments. Only for 
a few special cases {Di and D2 polynomials in m of degree smaller or equal than one and 
two respectively) these equations decouple. Otherwise, one has to rely on approximated 
methods to solve their dynamics. 

III. PARADIGMATIC MEAN-FIELD MODEL 

Let us exhibit our STD analysis for the paradigmatic system of A^ fully connected Ising 
spins (Curie- Weiss model), subject to a magnetic field H, ruled by the mean-field Hamilto- 
nian 

n = — —M^ - HM. (12) 

Since the Hamiltonian depends only on the total magnetization M, the master equation 
for this model can be written in closed form for M P, 112[. In the large N limit, when the 
magnetization per spin m = M/N can be taken as a continuous variable, an expansion of 
the master equation up to first order in the perturbative parameter e = 1/A^ leads for the 
Glauber dynamics to a FP equation ([7]) with [l2(] 

Di{m) = —m + tanh[m] — e(3Jmsech [m] , 

D2{m) = e(l — mtanh[r?7,']) , (13) 

where we have defined m' = P{Jm + if), with /3 = l/lksT). 

In the next sections we derive asymptotic solutions of the FPE with these coefficients, 
both close to the critical point [H = and T ^ Tc = J/ks) and to spinodal points for 
T < Tc. Analytical results are compared against Monte Carlo simulation ones using Glauber 
algorithm. Time was adimensionalized with the characteristic time to of the transition rate 



w = tQ^{l + exp{(3A'H))~^. The unit of time in theoretical expressions corresponds to one 
MC step in simulations. We also performed several checks using Metropolis algorithm. The 
outcomes were indistinguishable from the Glauber ones, except for a trivial time rescaling 
factor 2 close to the critical point, as expected TJ. 



IV. STD NEAR THE CRITICAL POINT 



In the vicinity of the critical point (at T c:^ Tc = J /ks = 1 and H = 0), the coefficients 
f lT3|) can be approximated for small m (i.e., /3J|r7i| << 1) respectively by 



L»2(m) = e ([1 - (1 - X)m^] + 0(m^)) , 



(14) 



where w(A, e) = A + e(l - A) and k{\, e) = (^ - e)(l - A)^, with A = 1 - TJT. 

Within the domain of validity of these approximations (1 — X)m'^ << 1 and therefore 
D2 — e. Concerning Di, its linear term dominates, that is. 



Di{m) ~ — w(A, e) m . 



(15) 
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FIG. 1: (Color online) Potential V{m), for different values of T in the vicinity of Tc, indicated in 
the figure. It was obtained (unless an arbitrary additive constant) from the integration of Di in 
Eqs. (|13p (black full lines), P^ (red dashed lines) and (fTCj) (gray daslied lines). Inset: zoom of 
the region close to m = 1. 



if 

|a;| >> Km^ . (16) 

This implies a parabolic approximation of the drift potential V{m) = — J Di{m)dm, whose 
shape is plotted in Fig. [T]for different values of T ~ Tc, found from the integration of Di in 
Eq. f lT3|) and of the linearized expression ([15]), for comparison. For u > 0, one has a confining 
quadratic potential, while for u < the parabolic potential is inverted, with an unstable 
point at m = 0. 

A. Ornstein-Ulhenbeck approximation 

Now, for linear Di and constant D2, the exact solution of Eq. ([7]) reads 

1 / Irri — mn exr)( —u t)]"^ s , _. 

where cr^(t) = e[l — exp{—2ujt)]/uj. This solution applies for a; > (Ornstein-Uhlenbeck 
(OU) process) as well as for u < 0, and is valid as long as the probability distribution 
remains strongly picked so that the inequality flT6l) holds for any value of m with non- 
negligible probability. 

Performing the average with Eq. ( TT71) gives 

(m) = moexp(— cut) . (18) 

Therefore, for u > (<)0, that is T/T^ > (<) 1 — 1/A^, the average magnetization decays 
(grows) exponentially, with characteristic time |cj|^^. Then, for time scales t ^ |i^|^^, it 
remains (m) ~ mo- Since in the large N limit w ~ A, then the magnetization scales as 
{m) = moF{Xt). This is consistent with Eq. ([2]), provided that ^ = and uz = 1, in 
agreement with the mean-field exponents u = 1/2 and z = 2. The same exponents are 
displayed by the Gaussian model 1|. For higher-order moments m'^"^ = {{m — (m))'^) with 
even n >2, one has 

^(n) ^ _±^^2eu-\l-exp[-2ut]p . (19) 



Then, for short times t << 1/\uj\, 



^W_[et]"/^ (20) 

Hence, m^'^^ ~ t/N, consistently with Eq. (j3]) (/3 = 1/2), provided that we choose d = A, the 
upper critical dimension. 



The characteristic time scale for STD behavior is then t ^ tstd with 

1 N 

Tstd 



If |A A^l ^ 1 we have tstd ~ Vl-^l ^ ^5 while for |A A^| <^ 1 we have tstd ~ N- 



(21) 
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FIG. 2: (Color online) First and second moments of the order parameter as a function of time 
t, for TTT-o = 0.01 and different values of T ~ Tc = 1. (a) Magnetization: black dot-dashed lines 
correspond to Eq. p^ and black dashed ones to Eq. (p2|) . (b) Second moment: black dot-dashed 
lines correspond to Eq. (J19p . Numerical simulations using Glauber dynamics were performed for 
A^ = 8 X 10^ (color symbols). 

Fig. 12] displays the comparison between numerical simulations and the approximate OU 
solutions Eqs. ([T8])-([l9]), for A^ = 8 x 10^ mo = 0.01 and different values of T ~ T^, such 
that |AA^| ^ 1. The OU approximation gives an excellent agreement for time scales up to 
t ~ Tstd [jstd ~ 100 for the present parameter values). Averages were taken over 1000 
independent MC runs. The main differences between the theoretical and numerical results 
appear for T < T^ and t > tstd, where finite-size effects shift the equilibrium value of both 
the average magnetization and its variance. 

Fig. [2] also shows the performance of Eq. ( 1221) . which reproduces the simulation results 
for longer times than Eq. (fTSl) . predicting the transient steady state. The lower saturation 
level observed in simulation outcomes for T < T^ is due to the presence of fluctuations 

8 



that drive some trajectories to the equihbrium state with negative magnetization, while the 
deterministic equation rules the stabilization at the level of the local minimum. Also notice 
that this discrepancy decreases as T departs from the critical value because of the consequent 
increase of the potential barrier height, which makes such events less probable. For T > T^, 
the system evolves quickly towards the vicinity of the equilibrium state and the saturation 
level of the second moment is very close to the value given by the (bimodal) steady state 
distribution P{m) oc exp(— l^(r?7,)/e). In any case, finite-size higher order corrections can be 
neglected as far as the STD behavior is concerned. 

B. Quartic approximation of the drift potential 



When (1161) does not apply, one can not discard the cubic contribution to Di. For such 
case we show in Appendix [A] that the inclusion of the cubic correction in the drift coefficient 
Eq. ([H]) leads for e < 1 to 

/ \ moe"'^* 

(m) = — = . (22) 

v/l + m2/t(l-e-2"*)/u; 

This solution is exact in the thermodynamic limit e — ?■ 0, as can be verified by direct 
integration of the deterministic version of Eq. flTTj) [7|, i.e., 

^ = fl.(W). (23) 

Notice that the expansion of Eq. ( l22l) up to first order in uiq reproduces Eq. ( TT8l) . The case 
u = {T = Tc), can also be drawn from Eq. fl22l) by taking the limit a; — )■ 0, yielding 

(-) = . To 2 , ■ (24) 

a/1 + 2mQ/tt 

In Appendix El we additionally show that finite-size corrections do not change the STD 
scaling of (m) . For the second moment we obtain 

m(2) ^ (m^) - (m)^ = 2et^i^^|^^^^^^ + 0{e^ eu) , (25) 

where z = urnQt. Notice that up to a typical time scale l/{2K,mQ), the approximation 
ijiC^) ~ 2et holds. For KuiQi >> 1, a crossover to a second linear (hence normal diffusive) 
regime but with a different diffusion constant is predicted, namely m^"^^ ~ ^i/'^, although it 
typically falls beyond the STD region. 
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FIG. 3: (Color online) First and second moments as a function of time t for rriQ = 1 and different 
values of T ~ Tc = 1. (a) Magnetization: dashed lines correspond to theoretical results given by 
Eq. (|22p . (b) Second moment: the full black line represents 2D2{mQ = l)t. Numerical simulations 
using Glauber dynamics were performed for N = 8 x 10^ (color symbols). 

C. Other initial conditions 



To investigate the scaling behavior for other initial conditions, we analyzed the STD 
behavior when mo = 1. As can be seen in the inset of Fig. [H the cubic approximation still 
holds close to m = 1. Hence, the thermodynamic-limit expression (12^ is expected to apply 
too, as verified in Fig. [3^. In comparison with the initial condition of Fig. |2l here trajectories 
get more trapped around the positive minimum, hence the agreement with deterministic 
Eq. ( 122|) is still better. For finite systems, the intensity of the fluctuations is state dependent 
following Eq. ( fT3l) . Therefore, the finite-size corrections derived by assuming D2 — e do not 
hold. However, for very short times one still expects rri^^-* ~ 2D2{'mo)t, according to Eq. (jQ]), 
as in fact verified in numerical simulations illustrated in Fig. [3l From Eq. ( I22l) we have 
that m(t) ~ t^^/^(l — At) for t <^ 1/1-^ 1 5 i^ agreement with Eq. ([5]). The excellent accord 
between Eq. (!22|) and numerical simulation outcomes displayed in Fig. [3] when |AA^| ^ 1 
confirms our previous assumptions. Numerical simulations for other values of A^ also verify 
the above scaling. For T > Tc the equilibrium (final steady state) values of both mean and 



10 



variance are quickly approached as in Fig|2l However, when T < Tc, we see from FigjSb 
that all the curves lie below the critical curve, at variance with the behavior observed when 
rrio ^ 1 (compare with Fig|2)D). This is because when rriQ = 1 almost all the trajectories 
get trapped in the positive minimum. Thus, the variance stabilizes in a value corresponding 
to the fluctuations in a single potential minimum. At long enough times, both minima in a 
finite size system get equally populated and therefore the equilibrium value of m^'^^ will be 
higher. However, the time scales needed to observe this effect fall outside the STD regime. 
On the contrary, when rrio <^ 1, a relatively large number of trajectories cross the barrier 
between minima and m'-^^ approaches the equilibrium value (which is larger than the steady 
one), even at very short times, as can be verified by comparing the numerical plateaux in 
FiglDD with the equilibrium value 

pi 2 -V(m) 

,^, /-.me E am 

< = I -_XS^, ■ (26) 

j_^ e <! dm 

V. STD NEAR THE SPINODAL 

When T < Tc the model has a line of first-order transitions aX H = and metastable 
stationary solutions for a range of values of H. Without loss of generality we will restrict 
hereafter to the metastable solutions with positive magnetization, that is, those analytic con- 
tinuations of the equilibrium magnetization from positive to negative values of H. Defining 
h = f3H, the metastable state exists as long as h > hsp, where the spinodal field is given by 

hsp = -pJmsp + - In 

2 1 — m^sp 



msP= ,11- jj 



le 



where m^sp is the magnetization at the spinodal point |6|. 

Suppose now that we start the system evolution from the completely ordered state m,Q = 1 
with T <T(. and h > hsp and let us define Am = -m — m^sp and A/i = h — hsp- Considering 
Am as an order parameter, numerical simulations using Metropolis dynamics [6J] showed 
that close enough to the spinodal point (|A/i| ^ 1) its moments obey the scaling form (jlj) 
with r = Ah/hsp- For temperatures far enough from Tc the spinodal magnetization m^sp 
is close to one and we can expand Di and D2 in powers of Ah and Am. Moreover, close to 



the spinodal we can neglect 12|] the finite-size correction of Di. Then, from Eqs. ( 1T3|) one 

11 
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FIG. 4: (Color online) Potential V{m), for different values of h in the vicinity of hgp, indicated 
in the figure for T = 4/9. It was computed (unless an arbitrary additive constant) from the 
integration of Di in Eq. ([13]) (red full lines) and (p7|) (black dashed lines). 

has at first order in Ah and second order in Am: 

Di(m) ~ — -— — 2mspA'mAh — (3Jmsp{AmY , 

D2{m) ~ e(^-2mspAm + {l3J-2){Anif 



(27) 



In Fig. Hlwe plot the shape of V{m) for different values of h in the vicinity of hsp-, obtained 
both from integration of Di in Eq. flT3|) and of the approximate quadratic polynomial ([2] 
for comparison. 

The moments of Am can be calculated by means of Eq. ([H]), namely 



((Am)") = Y,[DA + D2d,,fx^t^/k\ , 



(28) 



fc>0 



where we have defined x = m^ — msp. 

For n = 1 we can neglect in a first approximation the diffusion term, that is, at least for 
short times we can disregard finite-size effects. Then, from Eq. (1281) using Di = —A{x'^ + 
2Aax — a), with a = Ah/ {(3 J A) and A = pjm^sp one has (see Appendix [B|) 



= ^^ ^ + tanh( V7^t) _ ^^ 
^ ' ^ ' 1 + M tanh(^At) 

12 



(29) 
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FIG. 5: (Color online) Mean magnetization as a function of time t, for niQ = 1, T = 4/9 and 
different values of h. Black dashed lines correspond to the prediction given by Eqs. (j29p -(|3i p . The 
full black line corresponds to Eq. ()B4p . The inset is the same plot in linear-log scales. Numerical 
simulations were performed for N = 8 x 10^ (color symbols). 



where u = (x + Aa)/^/^ and -y = a + J^o? . For a < (hence 7 < 0), Eq. fl29|) becomes 



(Am) 



- Aa 



(30) 



1 + M tan(A/|7|v4t) 

Alternatively, Eqs. ( I29|) - (!30l) can be obtained by integrating Eq. (123|) . and are in good 
agreement with numerical simulations, as illustrated in Fig. ([5]). One observes the following 
asymptotic behaviors: 

(i) For \h\ < \hsp\ {a > 0), a constant level is reached. In fact, since the potential 
presents a local minimum, the plateau occurs at a level associated to that minimum. This is 
in accord with numerical simulations (Fig. [5]), notice that the local minimum of the potential 
is at m ~ 0.768, then Am = m — msp — 0.023, in agreement with the observed level. 

(ii) For \h\ > \hsp\: (a < 0), Eq. (I30ll yields a rapid decay towards zero attained at 
finite t. This is because the potential is tilted towards the absolute minimum (without local 
minimum) . 

In the limit a — ?■ 0, from Eq. (129|) it follows 

X 



(Am) 



(31) 



1 + Axt 

Hence, at the spinodal point one has (Am(t)) ~ t~^ for t ^ 1/Ax, consistently with Eq. i^ 
with (3 = 1/2 and uz = 1/2, in agreement with previous numerical results p]. This behavior 
corresponds to the relaxation towards the saddle point m = rrisp. While in an infinite system 
such point is an stationary state, finite-size fluctuations destabilize it, with the subsequent 



13 



exponential relaxation towards the equilibrium value Am ^ — 1 — hsp at longer times, as 
depicted in Fig. [51 Finite-size corrections to Eq. f l3T]) . that we compute for Ah = 0, can be 
obtained by including the diffusion term in Eq. fl28p . When Ah = 0, following Eq. fl27p . we 
have D2{x) ~ e(ax^ + bx + c), with a = PJ — 2, 6 = —2msp, c = 1//3J and Di = —Ax"^. 
In Appendix [B] we obtain Eq. (lB4p . furnishing (Am) corrected at first order in e, that for 
t ^ 1/Ax leads to 



(a™) ~ -L 



ecA 



2 



^~^^' + ^^'^''''^ 



(32) 
Hence, finite-size effects will become relevant only when t ^ t*, with 



1/3 /1n^AV3 



^--m-m^ 



in agreement with the scaling proposed in Ref. |6|: t* oc A^^/*^", with z = 2 and d^ = 6. 

Finally, let us consider the second moment. In Appendix |B]we obtain Eq. ( lB6p . that gives 
the e-correction to {{Amy). It allows to compute Am^'^^ = ((Am)^) — ((Am,))^, Eq. ( ]B7p . 
that at short times t ^ l/Ax leads to 

Am^^^ ~2e(ax^ + 6x + c)t ~2D(x)t, (34) 

in accord with Eq. ([9]). 

Meanwhile, for t ^ l/Ax, Eq. (1B6P behaves as 



{{AmY) 



1 



r\j 



{Atf 
Hence from Eqs. ( 132|) and (!35|) one gets 



ecA 



2 



,t3 + 0(et2^e2) 



(35) 



Am(2) ~ ^ . (36) 

Notice that in this regime, the pref actor of t given by Eq. ( 136|) is generically different 
from that obtained in the very short-time regime following Eq. fjMj) . Fig. ([6]) illustrates this 
cross-over for different values of ttlq and fixed temperature. The prefactor at small times, 
202(1710) varies with itlq (panel a), while at intermediate times 1/Ax ^ t < t* the prefactor 
becomes |ec = -^ independently of niQ, which is evident in the linear scale (panel b). 

In any case the behavior Am^'^' ~ et up to t ~ t* is consistent with the STD scaling 
hypothesis for the set of mean-field exponents zu = 2, P = 1/2, dc = Q and in agreement 



with numerical outcomes J6|. 
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FIG. 6: (Color online) Second moment of the order parameter as a function of time t, for T = 4/9, 
h = hsp and different values of mo (with tusp < mo < 1). Panels (a) and (b) display the same data 
in logarithmic and linear scales, respectively. Black full lines correspond to the prediction given 
by Eq. (|B6p . Symbols correspond to MC numerical simulations for A^ = 8 x 10^. The dash-dotted 
lines correspond to 2D2{msp)t (upper line) and 2/5 cet (lower line). 



VI. FINAL COMMENTS 

We studied the short-time dynamical behavior of finite-size mean-field models (infinite- 
range interactions) with non conserved order parameter dynamics. By solving the associated 
Fokker-Plank equation we obtained closed expressions for the first moments of the order 
parameter, in the vicinity of both the critical and spinodal points. This allowed us to confirm 
the STD scaling hypothesis in both situations, as well as to determine the dynamical ranges 
of its validity. In particular, we confirmed analytically its validity when the system starts 
from an ordered state. Moreover, we found that a diffusion-like scaling behavior of the second 
moment appears for any initial value of the order parameter, but the associated diffusion 
coefficient presents a crossover between two different values, for short and intermediate times 
within the STD regime. 

We found in general that the scaling behavior of the first moment is mainly determined 
by the shape of the potential V{m) = — j Di{m)dm and therefore by the equilibrium 
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generalized free energy f{m, T, H), which has the same extrenia structure as [T\ V{m). The 
scahng behavior of higher moments, on the other hand, has its origin on the Gaussian nature 
of finite-size fluctuations close to the singular points. Although our results were obtained for 
a particular model, it is worth to stress that the above facts are characteristic of mean-field 
systems, since they depend only on the shape of V{m) and on the proportionality D2 oc 1/A^. 
This makes the analysis quite general and independent of the particular mean-field model. 

Acknowledgments: The authors would like to thank T. S. Grigera and E. S. Loscar 
for sharing with us their simulation codes for Metropolis dynamics, as well as for useful 
discussions. This work was supported by CNPq and Faperj (Brazil), CONICET, Universidad 
Nacional de Cordoba, and ANPCyT/FONCyT (Argentina). 

Appendix A: Quartic potential approximation near the critical point 

To investigate the effect of including the cubic correction in the drift coefficient Eq. flT^ . 
we evaluate the particular setting of Eq. (EJ 

(m") = V[(-wmo - Kml)d^, + ed^,^,]''^— . (Al) 

fc>o '^■ 

In the limit |AA^| ^ 1, we can neglect in a first approximation the diffusion term and 
compute 



[y—UJiiLQ — i%iiLQ)Umo\ 
k>0 



(m) ^ ^[(-wmo - Kml)dmo 



k 



By iterating the operator k times and identifying the general form of the coefficients of t^, 
with the aid of symbolic manipulation programs, we obtain 

w«-E^^(p(-t)'E(:)Mr(-H-i)' 

k,j>0 J i=0 f- 






i>o J 



rriQe '^^ 



a/1 + mlK{l - e-2'^*)/w ' 
that coincides with the exact deterministic solution f l22|) . 



(A2) 
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Fluctuations can be neglected as long as k, w ~ Ole'^). However, while 3k remains of 
order one (except for extreme temperatures), typically a;~A + e<<l. Then, a finite-size 
correction can be included by keeping only the terms of order e and u in each coefficient of 
t^ in Eq. (lAip . This procedure yields the correction term, 

Ci(e) = -enmot^Yl ( i2k^ + 6k + 3){-z/2)'' 



-eKrriQt^ 



k>0 

3 + 42 + 2^2 

(1 + 2^)5/2 



where z = urnQt. Then, it results 



^""^ - (1 + 2^)1/2 +(i + 2z)3/2 y^^) 



^3 + 4:Z + 2z^ 



-enniot , + 0{e ,euj,u ) . 

(1 + 22;)^'^ 

Notice that the first two terms in the right-hand side come from the expansion of the 
deterministic Eq. ( ]A2|) up to first order in u. 

In particular, exactly at the critical point we have k, = 1/3 and A = (hence u = e). 
Therefore, as in the case of the OU approximation, one concludes that the magnetization 
remains m ~ rriQ up to a characteristic time tq ~ 1/e = A^. 

Similarly, for {m^), one obtains the correction 

C,{e) = 6t5^(A: + l)(A: + 2)(-2.)^ = -^|— , 

A;>0 ^ ^ 

(1 + 2z) k(1 + 2zY (1 + 2zY 
Since in the deterministic limit (m") = (m)"', then the first two terms in the right-hand side 
come from the expansion of the squared Eq. ( 1A2P up to first order in u. 

In the computation of the centered second moment, using Eqs. flA3P and flA4p . the purely 
deterministic terms cancels out to yield Eq. ( |25|) . 



leading to 
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Appendix B: Moment calculation near the spinodal 

If e = 0, from Eq. (^ using Di = -A(x'^ + 2Aax - a), a = Ah/((3JA) and A = (3Jmsp, 
the average magnetization is given by 

(Am) = 5^[-(x^ + 2Aax - a) d^]''x {Atf /k\ , (Bl) 

A;>0 

where x = tuq — msp. Completing squares and making the change of variables u = {x + 
Aa)/y/^ with 7 = a + A'^a'^ we obtain 

(Am) = y^^[(l - u^) dufu iVl^tf/kl - Aa . (B2) 

fc>0 

Considering the generating function for tangent, with the change of variable u = tanh^, 
one has [ij] 

^[(l-u2)aj"ur"/d = ^[9,]"tanhzr"/n! 

n>0 n>0 

= (u + tanhr)/(l + Mtanhr) , 

from where Eqs. (!29 l) -( !30l) follow. 

To include finite-size effects we have to consider the complete expression 

j-k 

(Am) = J2[Did. + D2d.^]'^ . (B3) 

When A/i = 0, from Eq. fl271) . we have D2{x) ~ e(ax^ + hx + c), with a = (3 J — 2, 
h = —2'msp, c = 1/(3 J and Di = —Ax^. The contributions of order e associated to each 
coefficient of the quadratic approximation of D2{x) are 

^^'^- 3A^^^ ^>^ y^ - 3A(i+i/)3' 

k>2 ^ ^' 



10(1 + y) 



where y = Axt. Summing the e-corrections Cia + Cib + Cic together with the deterministic 
one, given by Eq. ( 131]) . yields 



l + y Vl0yla;2 



,2 
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Likewise, we calculate 

((Am)2) = J2[Did. + D2d,.f^ . (B5) 



.x^t^ 



A;>0 

In this case, the contributions of order e are 



k\ 



^ 2ax sr^ I ^ + 2 1 / ,!, 2axy 



^^^ = -I^L(^ + l)(^ + 2)(3fc + l)(-^) = QA{l + yY ' 

C. ^ c.(2 -, 1 X:(. + !)(. + 2)(2. + 1)(-.)-^) ^ ct(10 + 10.H-10.J^^+V^,.) ^ 

Summing up the corrections C2a + C'2b + ^^20 together with the deterministic term (given 
by the squared Eq. flM]) ). yields 

2 

((Am)2) = ^-^ + (^^(10 + 10y + 101/2 + 5y3 + 2/4) 

^/.o , -... , ...2 , „,3N , o_.A ey 



+ -(12 + 6y + V + r) + 2ax 1 ^^^ . (B6) 

Finally, the second moment is obtained through Am'-^-* = ((Am)^) — ((Arri))^. The purely 
deterministic terms cancel out and at first order in e it remains 

+ 156x(2 + y)(2 + 2y + y2) + 

+ 12c(5 + lOy + lOy^ + hy^ + y^)] + 0{t^) . (B7) 
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